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We study dynamical mass generation in QED in (2+1) dimensions using Hamiltonian lattice methods. We use 
staggered fermions, and perform simulations with explicit dynamical fermions in the chiral limit. We demonstrate 
that a recently developed method to reduce the fermion sign problem can successfully be applied to this problem. 
Our results are in agreement with both the strong coupling expansion and with Euclidean lattice simulations. 



1. INTRODUCTION 

QED in (2+1) dimensions exhibits both dy- 
namical mass generation @,|,|,U| and confine- 
ment This makes it an interesting theory 
to study these nonperturbative phenomena. Fur- 
thermore, these phenomena in QED3 are relevant 
for applications in condensed matter physics, in 
particular in connection with high-T c supercon- 
ductivity 0). Here, we focus on the question of 
chiral symmetry breaking in QED3. 

We use staggered fermions, and perform sim- 
ulations with explicit dynamical fermions in the 
chiral limit. Simulations with explicit fermions 
could improve our understanding of the micro- 
scopic behavior of (chiral) fermions in lattice 
simulations. We apply a recently developed 
method M , called the "zone method" , to reduce 
the fermion sign problem in our simulation. Pre- 
liminary results for the chiral condensate in QED3 
are in agreement with other methods. 

2. QED IN (2+1) DIMENSIONS 

QED3 is a superrenormalizable theory, with a 
dimensionful coupling: e 2 has dimensions of mass. 
This dimensionful parameter plays a role similar 
to Aqcd in QCD. In the chiral limit, it sets the 
energy scale. 

We use 4-component spinors, such that the 
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fermion mass term is even under parity. With 
one massless fermion, the Hamiltonian exhibits a 
global U(2) "chiral" symmetry. A fermion mass 
term breaks this symmetry to a U(l) X U(l) sym- 
metry The question is: is this chiral symmetry 
broken dynamically? The order parameter for 
this symmetry breaking is the chiral condensate. 

There are extensive studies using the Dyson- 
Schwinger equations suggesting that there is 
dynamical chiral symmetry breaking in QED3 if 
the number of fermion flavors is smaller than 
some critical number iV c ~ 3.3. However, the 
scale of this symmetry breaking (i.e. the magni- 
tude of the condensate) is rather small. 

2.1. Hamiltonian Lattice Simulations 

Following the notation of Ref . j| , we consider 
the lattice Hamiltonian in (2+1) dimensions 
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where r)x(f) = (-l) r2+1 , r] 2 (r) = 1, y = 1/g 2 , 
and Up is the usual plaquette operator. We use a 
dimensionless mass parameter \i = 2m/e 2 and a 
dimensionless coupling constant g 2 — e 2 a where 
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a is the lattice spacing. The continuum limit cor- 
responds to y — > oo or g 2 — > 0. 

We use staggered fermions, corresponding to 
one 4-component fermion flavor in the continuum 
limit [§J§). For /i = 0, the lattice breaks the 
"crural" U(2) symmetry to a discrete symmetry 
generated by a shift of one lattice spacing. A 
nonzero mass term breaks this discrete symme- 
try. To study chiral symmetry breaking on the 
lattice, we calculate the lattice condensate {'tpip) 
in the chiral limit as function of the lattice cou- 
pling y. It is related to the continuum condensate 

y 2 (^) utticc = ^(V^) continuum , (5) 

in the limit y — > oo. 

We use a finite spatial grid, Ni x N 2l and cal- 
culate observables 

, K = (^finalle^Ce^ | ^initial) 

{ ' ~ (*fi„al|e- fft ^initial) ' ( ' 

by dividing the Euclidean time variable t in N t 
slices, and inserting a complete set of states be- 
tween the time slices. We choose a tensor product 
basis of gauge and fermion states. We use a stan- 
dard Metropolis algorithm to update the gauge 
field configurations. 

For the fermions we use the worldline for- 
malism in combination with the loop algo- 
rithm to update the fermion configurations. 
As initial/final states we choose the strong cou- 
pling (y — > 0) ground state. The initial and final 
fermion states are identical, up to permutations 
of fermions. Numerically, these permutations of 
course give rise to the fermion sign problem. With 
the worldline formalism we can keep track of the 
permutations, and we use the newly developed 
zone method § to reduce the sign problem. 

2.2. Zone method 

The basic idea || is to introduce an n\ x n 2 
spatial sub- lattice (with m < Ni), which we call 
a zone, for the fermion simulations, see Fig. |l|. 
Next, we allow for permutations of initial and fi- 
nal state fermions inside this zone only. We do 
not allow for any configuration with fermion per- 
mutations outside the m x n 2 spatial sub-lattice. 
Note that at intermediate time slices the fermions 



Figure 1. Illustration of the zone method: an 
6x8 spatial lattice with a 3 x 4 zone. The size 
of this zone is 17, as characterized by the number 
of links inside the sub-lattice. 



are still allowed to wander through the entire 
Ni x N 2 lattice. In order to further reduce the 
fermion sign problem, we also use a larger value of 
the coupling y' and a larger value of the fermion 
mass m! when the fermions are outside this zone. 

It turns out that observables scale linearly in 
the zone size, provided that this zone size is 
larger than the characteristic "fermion wandering 
length" H . For any finite values of the coupling y 
and of the time variable t, this wandering length 
is finite, even for massless fermions. Thus one 
can extrapolate the results from relatively small 
zones to the entire lattice. 

3. NUMERICAL RESULTS 

There are different ways to define the size of a 
zone: by number of lattice points or by the num- 
ber of links inside the zone. For large lattices it 
does not matter which one one chooses. However, 
for the relatively small lattices we have used sofar, 
it turns out that the best way to characterize the 
zone size is the number of links inside the zone. 

In Fig. U we show that, within error bars, the 
chiral condensate indeed behaves linearly in the 
number of links inside the fermion zone, even for 
the smallest zones (at least for these parameter 
values). Thus we can use a linear extrapolation to 
obtain the lattice condensate for the entire lattice, 
using three or four different zone sizes. 

In Fig. U we show our results for a range of val- 
ues of the coupling y. These results are obtained 
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y'=0.15, \i 
•y' = 0.15, \i' = \i 
y'=y, = 1.0 
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Figure 2. Numerical results for different zone 
sizes within a 6 x 6 lattice, for fixed values of 
y = 0.5, (i = 0, t = 1.5, and N t = 10. The to- 
tal number of links in the lattice is 72, where the 
three linear fits meet (within error bars). 



on lattices up to 8 x 8 using several different zone 
sizes and for finite values of t. The error bars 
are statistical error bars only. We expect the fi- 
nite size effects due to the spatial lattice to be 
smaller than these error bars; however, prelim- 
inary results indicate the t — oo extrapolation 
could lead to a reduction by about 10%. We are 
currently working on estimating the finite size ef- 
fects in more detail [|l2[ . 

For y < 1 these results agree quite well with the 
strong coupling expansion [|J , whereas the strong 
decrease with increasing y for 1 < y < 2 is in 
good agreement with Euclidean Monte Carlo sim- 
ulations ||. It is unclear yet whether or not the 
condensate obtained this way is small or exactly 
zero in the continuum limit, y — > oo. 
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